Quantum noise and self-sustained radiation of T'T-symmetric systems 
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The observation that PT-symmetric Hamiltonians can have real-valued energy levels even if 
they are non-Hermitian has triggered intense activities, with experiments in particular focusing on 
optical systems, where Hermiticity can be broken by absorption and amplification. For classical 
waves, absorption and amplification are related by time-reversal symmetry. This work shows that 
microreversibility-breaking quantum noise turns PT-symmetric systems into self-sustained sources 
of radiation, which distinguishes them from ordinary, Hermitian quantum systems. 

PACS numbers: 42.50.Lc, 03.65.-w, 03.65.Nk 
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A common factor in quantum systems with a non- 
Hermitian Hamiltonian is the non-conservation of par- 
ticle number, either because the system is leaky, or be- 
cause there is loss or gain in an absorbing or amplifying 
medium. Ignoring nonlinear effects such as the feedback 
in a laser, such systems ordinarily do not possess sta- 
tionary states; instead, they only support decaying qua- 
sibound states with complex energy, where the imagi- 
nary part Ini-E = — l/2r (setting h = 1) accounts for 
particle loss with decay rate 1/r (particle gain is asso- 
ciated to a negative decay rate). A notable exception 
are non-Hermitian systems with loss and gain combined 
such that they are invariant under joint parity (P) and 
time-reversal (T) symmetry If there is no leakage, 
these PT-symmetric systems generically possess a set of 
real-valued energy levels, as well as complex-conjugate 
pairs of complex energy levels. Systems with entirely 
real spectrum define a consistent unitary extension of 
quantum mechanics 0, Hj]. This observation has led to 
intense research efforts delivering a new theoretical per- 
spective on systems as varied as quantum field theories 
and complex crystals [H , while experimental realizations 
in particular focus on optical systems where Hermiticity 
can be violated by absorption and amplification [5[. 

For classical waves, amplification and absorption are 
strictly related by time reversal. The existence of sta- 
tionary states with real energy can therefore be seen as a 
consequence of the balance of amplification and absorp- 
tion in parity-related regions of a 'PT-symmetric system. 
At the heart of absorption and amplification, however, 
are noisy microscopic quantum processes (such as spon- 
taneous and stimulated emission events, and stimulated 
absorption events) which effectively break time-reversal 
symmetry Q. The objective of this Letter is to show 
that the effects of this quantum noise distinguish PT- 
symmetric systems from Hermitian quantum systems, 
and indeed suggest an alternative interpretation of the 
physics behind non-Hermitian PT symmetry: (i) Ac- 
counting for quantum noise, PT-symmetric systems with 
stationary states are self-sustained sources of radiation, 
fed by the pumping in the amplifying parts of the system, 
(ii) That the energy of these states is real means that the 
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FIG. 1: Illustration of the scattering input-output approach 
to non-Hermitian PT-symmetric systems, defining the scat- 
tering input-output operators S^'r"* and internal bosonic 

noise operators bL,n in different parts of the system. The op- 
erators aQ lght ' lc 4 and a" 1 ^ "' feature in the intermediate steps 
of the calculation (for details see the text). Semitranspar- 
ent mirrors with transmission probability F are introduced to 
study the limit F — > of a closed (non-leaky) system. (In 
the context of present experiments with optical fibers 0] , this 
sketch relates to the transverse confinement.) 



system is stabilized at the lasing threshold, (iii) When 
the system is leaky, the emitted radiation breaks parity 
symmetry (i.e., the emission pattern is asymmetric), (iv) 
In the limit of a non-leaking system, the emitted radia- 
tion intensity approaches a constant value, and provides 
a direct measure of the non-Hermiticity of the system. 
The internal energy density of radiation then diverges, 
which entails a practical limitation for the implementa- 
tion of PT symmetry in closed systems. 

These conclusions are obtained by employing the 
quantum-optical input-output formalism [7[ in its scat- 
tering formulation [8rllol|. The scattering approach also 
provides insight into PT symmetry for classical waves 
[Til [T3| . which defines the starting point of this Letter. 

Scattering approach to non-Hermitian VT -symmetric 
systems. — Probing the internal dynamics of an optical 
system by external radiation naturally leads to the scat- 
tering scenario depicted in Fig. [TJ The relation a out — 
Sa ln between incoming and outgoing wave amplitudes is 
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provided by the scattering matrix 
S(E) = 



r t' 
t r' 



(1) 



which contains blocks describing reflection (r, r') and 
transmission (t, t') when probed from the left or right, re- 
spectively. Each block consists of an TV x iV-dimensional 
matrix, where N is the number of modes at each en- 
trance. The scattering matrix is generally energy depen- 
dent (which arises from the wavcnumbcr dependence of 
constructive and destructive interference), and its poles 
determine the energies of quasibound states in the system 
0. 

In general, the scattering matrix fulfills the following 
two reciprocity relations: The Onsager relation 



S( 1 ,-B,E) = S T ( 1 ,B,E), 



and the relation 



S(-j,B,E) = [St ( 7j B , E*)Y 



(2) 



(3) 



of classical microreversibility. Here, 7 and B characterize 
two possible sources of broken time-reversal symmetry: 
absorption or amplification (7 > or 7 < 0), which con- 
tribute an imaginary symmetric (non-Hermitian) term 
to the Hamiltonian H, and magneto-optical effects (B), 
which contribute an imaginary antisymmetric (but still 
Hermitian) term. 

Conventional time-reversal TH = H* transforms so- 
lutions according to 1~iJj — tp* , which interchanges in- 
coming and outgoing states, and therefore transforms the 
scattering matrix according to 

TS(j,B,E) = [S*(^,B, E)]- 1 = S(- 7 ,-B,E*). (4) 

Assuming that energy is real, a system has T symme- 
try, TS = S hence S* = S' 1 , if S(-/,B) = S(-j,-B), 
which requires 7 = B = 13|. Parity VH{x) = H(—x) 
transforms solutions according to Vij)(x) = ip(—x), which 
exchanges the left and right leads and yields 



PS(~,,B,E) = a x S(j,B,E)a x , 



(5) 



where a x is a Pauli matrix operating on the blocks in 
Eq. ([T]). The VT operation on the scattering matrix is 
therefore given by 

VTS(j,B,E) = o- x [S*(7,fl,£7)]~V B (6) 
= a x S(-j,-B,E*)a x . (7) 

For Hermitian systems, VT symmetry implies S = 
a x S T a x [l5| . For non-Hermitian systems, VT symme- 
try implies the additional condition V~f = —7 [7(2) = 
— 7(— x)], i.e., there is a balance of absorption and am- 
plification in parity-related regions. 



Let us now explore from the scattering perspective how 
real-energy bound states appear in 'PT-symmetric sys- 
tems. As shown in Fig. [TJ such systems can be con- 
structed by joining two regions, where the left region, 

n t[ 

h r[ 

to the right region, S2 = VTS±, which using standard 
block-inversion formulas can be written as 



with scattering matrix Si = 



, is 'P'T'-symmetric 



/ 



S 2 



(ri-tirf 1 ^)' 



1 



\ 



(txrrX-ri)* 



(n - t'^-Hi)* I 

(8) 



Bound states can be studied by closing the system off 
by mirrors with small transmission probability r « 1, 
described by a scattering matrix 



Sv 



v 7 ! - T iVT 



(9) 



The mirrors can be included by wave matching at their 
inward-facing faces, which amounts to an algebraic elim- 
ination of the amplitudes a^' 2 ont in Fig. [T] 

The scattering matrix of the left half of the system 
then reads 



Sl 



\ 



1 + 



V i + r lv / T Tr r 1 + T 



/ 

(10) 

while the scattering matrix Sr = VTSl of the right half 
again follows from symmetry. These scattering matrices 
relate amplitudes of ingoing and outgoing modes (defined 
in Fig. [1]) according to 



right 



,lcft 



Joft 



'lit 



(11) 

The scattering matrix of the composed system is ob- 
tained by algebraically eliminating the amplitudes ag eft 
and aQ lght at the interface between both regions. For 
r — > 0, these amplitudes become singular when 



det Im (r' L ) = det 



Im 



tit[_ 

n 



1 



0, (12) 



which (due to the general connection of scattering poles 
and quasibound states) is the quantization condition of 
the closed system. Determinantal quantization condi- 
tions of this kind have been introduced for general sys- 
tems in Ref. 14|; they also form the basis of exact nu- 
merical techniques as reviewed in (l6| . The 'PT-specific 
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version (IT21) of the quantization condition requires that 
the N real column vectors of Im(r^) be linearly depen- 
dent, which generically can be achieved by varying a 
single real parameter such as E (identifying this as a 
problem of codimension one). Therefore, we recover that 
closed 'PT-symmetric systems typically possess a num- 
ber of bound states with real energy, even though the 
Hamiltonian is not Hcrmitian. Condition (|12[) also ad- 
mits complex eigenvalues, which then appear in complex- 
conjugated pairs. 

Quantum noise. — The scattering approach can be ex- 
tended to include quantum noise by passing from wave 
amplitudes a" 1 , a out to bosonic annihilation operators 
a m , a out , respectively. This defines the scattering vari- 
ant of the input-output formalism , which has been 
used to describe systems that are exclusively absorbing 
or amplifying. To adapt the approach to PT-symmetric 
systems, where both effects are combined, we formally 
separate the absorbing regions from the amplifying re- 
gions, and then join them together similar to the de- 
scription of classical waves, given above. 

For dcfinitcncss let us assume that the left half of the 
system is purely absorbing. For this part, the input- 
output scattering relations then take the form 



- right 



= Sr 



r.lefi 



Qhbh, 



(13) 



which connects the ingoing and outgoing modes to 
bosonic operators and b R representing quantum fluc- 
tuations in the left and right part of the medium. These 
operators can be defined microscopically in the frame- 
work of system-and-bath approaches, or phenomenologi- 
cally as operator- valued Langevin forces 0, H[ . The prop- 
erties of these operators can be determined from the con- 
dition that both a ln and a out have to satisfy standard 
canonical commutation relations. This dictates that the 
coupling matrix Ql satisfies the fluctuation-dissipation 
theorem Q 



QlQI = i - s L sl 



(14) 



In the right half of the system, where the medium is 
amplifying, we have 



Aleft \ / -right \ 

Sut )=Sr (% n )+Q R b R , (15) 



where the commutation relations now dictate coupling to 
creation operators, with 



QrQr — SrSr !• 



(16) 



By assumption, the operators b\ and b R commute with 



a m ; however, according to Eqs. (|T3|) and (fT5|) they do not 
commute with a out , which is a manifestation of broken 
micro-reversibility in quantum optics. 



We can now describe the full PT-symmetric system by 
algebraically eliminating the interface operators fig and 
clq s . In the absence of incoming radiation, the intensity 
emitted to the left and right is then given by 

I L (E) = ±-(a°^ar), lR(E) = ±(aT^T), (17) 



which can be evaluated assuming 

(b{b L )=0, (b R b R } = 



(18) 



(ground-state population in the absorbing regions and to- 
tal population inversion in the amplifying regions; these 
conditions minimize the quantum noise). 

Let us first consider the case of a single-mode resonator 
(N = 1) with purely ballistic internal dynamics and ab- 
sorption in the left region [ItJ , with scattering matrices 



Si = 



o h 
h o 



S 2 = 



l/t\ 

i/t; o 



(19) 



where \t\\ < 1. Including the mirrors, the total scattering 
matrix is UM 



( VT^T{tf -t\) 



S 



M 2 r 



\ 



t((i 



(i-r)-t; 



M 2 r 



t\{i -r) -tf 



vT=r(fi 



*2 



(20) 



V t\{l-T)-tf tf(l-T)-tf I 



and the quantization condition (|T2|) for the closed res- 
onator takes the form Imtf = (indeed, the scatter- 
ing matrix elements diverge under this condition and 
r — > 0). Following the quantum-optical procedure de- 
scribed above [and using, in particular, Eqs. (fT4")l , p^|) . 
and qi8|]. we find that this resonator emits radiation of 
intensity 



Il(E) 



r(|*i 



-l l-r 



M 2 ) 



2 7 r|(i 1 /tt) 2 -l + r|2 

r(i-M 2 ) (i-r + i^r 2 ) 



27r|(tx/^) 2 -l + r|2 
Since \ti\ < 1 this gives Ir > II, the difference being 



(21) 
(22) 



AI(E) = I R (E) - I L (E) = 



r 2 (\h 



\ti\Y 



2 7 r|(t 1 /^) 2 -l + r| 



(23) 



Therefore, the emission from the right exit, close to the 
amplifying region of the medium, is larger than the emis- 
sion from the left exit, close to the absorbing region of 
the medium. The overall output intensity to both sides 
can be written as 



I(E) = I L (E) + I R (E) = 



r(2-r)(|^ 1 |- 2 -|t 1 | 2 ) 
27r|(*iAi) 2 -i + r| 2 ' 



(24) 
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Close to quantization in the closed system [T <C 1, 
E w Eq, where Eq fulfills the quantization condition 
Im^(_E ) = 0], the emission pattern becomes symmetric 
(as a consequence of VT invariance) and approaches a 
Lorentzian of the form 



Il(E)=I r (E) 
Here to = ti(Eo), while 



r(M- 2 -M 2 ) 

2Tr\2ir(E - E ) +T\ 



21m- 



dtx 
dE 



2L/c 



(25) 



(26) 



E—Eq 



is the transmission delay time of propagation between 
the two mirrors (with L the length of each region and c 
the speed of light) . The full width at half maximum is 
given by AE = T/t. While this width shrinks to zero as 
the system is closed off, remarkably the total intensity 



1(E) dE 



\to\- 2 -\to? 
2t 



(27) 



remains finite, and can be interpreted as a direct mea- 
sure of the degree of non-Hermiticity of the system (for 
ballistic transport, Hermiticity implies |<o| = !)■ 

In the case of a single-mode resonator with backscat- 
tering (where n and are finite), compact expressions 
can still be obtained as long as the leakage remains 
small (r <C 1), implying according to Eq. (fTU|) that 
\tl + 1|, | til, \ t' L \ <C 1. The emission pattern is then still 
symmetric, with intensity 



Il(E)=I r (E) 



1 



(1 



r-i |2M+> 12 



2tt |2(Imr^) 



(28) 



Linearization around the quantization condition again re- 
veals a Lorentzian line shape, with line width 



AE = Re 



d[(lmr' L )/t L t' L 
dE 



(29) 



Accounting for the scaling (ITU1) of scattering coefficients 
with L, the total intensity /tot oc (1— l^l 2 ) again remains 
finite as L — > 0. In the Hermitian case, this limit would 
imply \r' L \ — 1, so that the intensity vanishes. There- 
fore, the emitted radiation is still a direct measure of the 
degree of non-Hermiticity. 

Following the general formalism described above, the 
observations for one-dimensional scattering can be ex- 
tended to the general case of 'PT-symmetric systems with 
many modes, for which compact expressions are no longer 
available. The emitted intensity generally remains finite 
even in the limit of a closed system. Because the expecta- 
tion values (4 cftt a[, cft ) cx r-\ (a^ ightt ^ ight ) « F- 1 of the 
internal operators formally diverge in this limit, this is ac- 
companied by a diverging internal energy density, which 
can be interpreted as the source of this radiation. De- 
viations from perfect VT symmetry can be incorporated 



into the scattering formalism by taking Sr ^ VTSl, 
which results in an asymmetric emission pattern even in 
the limit of non-leaky mirrors. 

In standard laser theory [l9| , an active optical medium 
is below threshold (and stable) as long as all fundamen- 
tal modes are damped, which is characterized by complex 
energies with a negative imaginary part. With increasing 
pumping, the loss of the modes within the amplification 
window is counteracted by gain, and as soon as one of 
the energies acquires a positive imaginary part the sys- 
tem becomes instable: the intensity increases exponen- 
tially until saturation sets in. Therefore, the radiation 
features described above are characteristic of a laser sta- 
bilized at threshold, the threshold condition being that 
the energy of the lasing mode is real. The system then 
is marginally stable, and the internal radiation energy 
accumulates due to the quantum fluctuations. In prac- 
tice, this leads to saturation in the amplifying parts of 
the system and therefore identifies an obstacle for the 
implementation of strict VT symmetry in closed optical 
systems. When the system is slightly opened up, the 
threshold condition is no longer met and the internal en- 
ergy density is finite, while the emitted intensity remains 
a direct measure of the non-Hermiticity of the system. 

Concluding remarks. — In summary, the present Letter 
demonstrates that accounting for quantum noise, opti- 
cal realizations of non-Hermitian PT-symmetric systems 
emit radiation of an intensity that provides a direct mea- 
sure of non-Hermiticity. 

This has both practical as well as fundamental impli- 
cations. The fundamental consequences arise because 
non-Hermitian PT-symmetric systems with an entirely 
real spectrum define a consistent unitary theory of quan- 
tum mechanics 0, Q , formalized by the concept of quasi- 
Hermiticity, which introduces a new scalar product based 
on a generalized conjugation operation C that satisfies 
C 2 = 1, [C,H] = [C,VT] = 0. The self-sustained radia- 
tion identified here shows that accounting for quantum 
noise, non-Hermitian PT-symmetric systems are physi- 
cally distinct from ordinary Hermitian quantum systems: 
the canonical commutation relations for the input and 
output operators are only invariant under unitary trans- 
formations, which constraints the possibility to introduce 
alternative scalar products. This enforces classical dis- 



tinctions based on the transparency of such systems [12 1 . 

From a practical perspective, the self-sustained radi- 
ation can be used as an indicator of successfully imple- 
mented non-Hermitian VT symmetry in leaky systems, 
while the accompanying marginal instability and diverg- 
ing internal energy density signifies a practical obstacle 
for its implementation in the limit of no leakage. Fur- 
thermore, these findings identify a hitherto unexplored 
arena to study quantum fluctuations in active systems 
close to the lasing threshold. 

Non-optical realizations of PT-symmetric systems of- 
fer a wide range of connections to quantum field theories 
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[J], and while the results in the present paper are not 
directly transferable in detail, they suggest that quan- 
tum noise should provide fundamental insight into such 
systems, as well. 
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